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Abstract. - We study the high-velocity-projectile limit of the polarizational contribution 
to the in-plane stopping power in a strongly coupled two-dimensional electron liquid. The 
dielectric formalism based on the method of moments is employed. The frequency moments 
of the loss function are calculated using the model Hamiltonian including the two-dimensional 
Coulomb interaction potential proportional to the inverse power of k. We prove that the leading 
term of the high-velocity asymptote, like in the random-phase approximation, is not affected 
by correlations. 



In the last two decades, there has been a great interest in the investigation of static and 
dynamic properties of two-dimensional charged particle systems. Two-dimensional electron 
liquids trapped above the free surface of liquid helium [1] constitutes an example of this kind 
of systems within the classical regime. 

Moreover, recent advances in modern semiconductor nanotechnology have made it possible 
to fabricate quasi-two-dimensional quantum well structures of electron layers conforming a 
degenerate strongly correlated liquid phase. This is the case, for instance, of the systems 
consisting of electrons confined in the vicinity of a junction between a semiconductor and 
insulators [2], e.g., in a metal-oxide-semiconductor field-effect transistor, MOSFET, structure. 
Here, quasi-two-dimensional means that electrons have quantized energy levels along one 
dimension, whereas they are free to move within the two-dimensional junction. 

The system in question is a two-dimensional negatively charged one component plasma 
(OCP) immersed in a uniform and rigid neutralizing background of positive charges. We 
presume that electrons interact via a two-dimensional Coulomb-like potential, inversely pro- 
portional to the distance. 

There are two dimensionless parameters, which usually characterize these systems: the 
first one is used to quantify its degeneracy: D = PE-p, where /3 is the inverse temperature 
in energy units and Ep = h^k^/2m is the Fermi energy, fcp = V^Tm being the Fermi wave- 
vector, m the electron mass and n stands for the two-dimensional electron density; the second 
dimensionless parameter is the coupling parameter, F = /3e^/a, where a = (irn) is the 
2D Wigner-Seitz radius and e the electron charge. Whereas the parameter F characterizes 
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classical Coulomb systems, for quantal systems we might replace the temperature, f3~^, by 
Ep. In this case F equals the Brucckncr parameter, rg = a/ae, an being the Bohr radius. 

Two-dimensional electron liquids arc known to crystallize at about F » 137 ± 15 in the 
classical case [3], and at rg ~ 37±5 for the zero-temperature system [4], although this value can 
be affected by the presence of disorder [5]. For Coulomb systems, the strongly coupled regime 
is characterized by a value of the coupling parameter much higher than unity. Nowadays, this 
kind of experimental physical conditions have become usual. In this situation, interparticle 
correlational effects become crucial and, consequently, to describe such systems theoretically, 
we need to go beyond the random-phase approximation (RPA). 

In general, interaction of charged particles with condensed matter can be studied by means 
of the system stopping power. This function might be helpful in order to derive general 
diagnostic tools for this kind of physical systems. The stopping power accounts for the energy 
lost by an external positively charged projectile (which is assumed to move along a classical 
trajectory) as it passes through matter. There are several ways of dealing with this problem. 
Within the linear dielectric approach [6, 7] one assumes the incident projectile to have a 
constant velocity, and its energy loss is related to the strength of the dipole created between 
the ion and the center of charge of the induced charge density. This physical picture becomes 
more accurate for high- velocity and heavy projectiles. On the other hand, under the scattering 
picture [8-10] one relates the stopping power to the transfer of momentum produced by 
electron collisions with an effective potential created by the moving ion. Further works have 
dealt also with nonlinear screening effects by means of the quadratic response approach within 
the RPA [11] and the employment of the density functional theory [12]. 

In the sequel we will focuse on the part of the energy losses due to linear polarization of 
the medium. It was Lindhard [6] who derived a general expression relating the imaginary 
part of the inverse dielectric function to the polarizational losses valid in three-dimensional 
systems. Thus, in this context, the derivation of a consistent dielectric formalism has played 
a key role. 

More recently, Bret and Deutsch [7] obtained a similar expression for the in-plane polariza- 
tional stopping power in a two-dimensional Coulomb system. For a point-like charged particle 
this expression reads as 
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Moreover, these authors derived a dielectric function within the RPA and obtained asymp- 
totic expressions for the polarizational stopping power for high- and low-velocity projectiles. 
In the former case, in contrast to the well-known 3D asymptote [13 15], it was shown in 
Ref. [7] that the leading term of the stopping power is inversely proportional to the projectile 
velocity. This result has been later confirmed by means of different approaches [8-10, 16, 17]. 
The influence of electronic correlations in a zero-temperature electron gas on the low- velocity 
asymptote was also considered by Wang and Ma [16] who used a parametrized expression for 
the static local-field correction factor. However, these authors explicitly neglected the effect of 
electronic correlations on the high- velocity limiting stopping power, an assumption which was 
entirely based on a physical picture. Our goal here is to provide rigorous argument supporting 
this physically sound assumption. 

More precisely, our dielectric formalism, which is based on the method of moments, permits 
us to account for correlational contributions and handle the problem analitically in order to 
show that the first term of the asymptotic expansion of the two-dimensional stopping power 
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is unaffected by electron correlational effects. 

Dielectric formalism. - The (truncated) method of moments provides a suitable mathe- 
matical tool in order to reconstruct response functions or susceptibilities of a non-perturbative 
physical system from its known frequency (power) moments or sum rules [18]. In particular, 
for the inverse dielectric function we refer to the frequency moments of its imaginary part (we 
assume the system to be isotropic), 

1 1"°° 

Cy^ / uo^-^luve'^ {k,uj)Auj, v = Q,l,.... (2) 

Tl" J-oo 

Here, odd moments (;/ = 2n + 1, n = 0, 1, . . .) vanish due to the parity of the loss function, 
— Im£~^ {k,u>) /uj. In addition, only the moments Cq, C2, and C4 converge, in general. These 
are related to the static properties of the electron system [18, 19]: 

Co(fc) = l-£-i(fc,0), (3) 

C2(fc) = ^fc^^L, (4) 
W2D = i^2Ti (k) being the 2D plasma frequency, 

C4{k)=u;l^{l + K{k) + L{k)), (5) 

where 

Kik)^l^k+-^e (6) 

is the kinetic contribution to the fourth sum rule, -Bkin being the average kinetic energy per 
electron, which, for example, will coincide with /3~^ for a classical system. Furthermore, the 
correlation contribution reads as 

^(^)-^ E ^^[^(k + q)-^(<7)], (7) 

which is applicable for the two-dimensional Coulomb system. It must be stressed that, whereas 
the sum rules C2 and C4 are obtained by the usual derivation procedure from the Kubo 
formula, the moment Co stems from the Riesz-Hergiotz transformation [20, 21] for a non- 
magnetized plasma [18]. 

Further, owing to the Kramcrs-Kronig relations and the Nevanlinna formula [20,21], the 
expression 



-W, ^ ^ , , ^2D (fc) [Z + Q (k, Z)] , . 



establishes a one-to-one correspondence between the set of all possible dielectric (complex) 
functions satisfying the sum rules (l3][5|), and an arbitrary parameter function q{k,z). The 
latter must be analytic in the upper half-plane Imz > and possess a positive imaginary part 
there, and, in addition, it must fulfill the condition: lim^^oo 1 [k, z) / z = for 9 < arg(z) < 
TT-e {0<9 < 7r/2). 

The system parameters appearing in ([8|) depend only on the frequency moments: 

ivf (k) = C2/C0 = ioi^ (fc) (1 - e-' {k,0)y' , (9) 
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ul (fc) = C4/C2 = c^Id (A:) {l + K{k) + L (fc)) . (10) 

The function q(k,z) is phenomenologically related to the damping of the plasmon mode 
of the system. Actually, if we put q {k, z) — iO"'', we obtain for the loss function: 



Ime ^ (fc, Lo) 



^2D 

2ujI 



5 [u + UJ2) + 



^-Co5{lo) + '^8{uj-lo2) 



(11) 



which resembles the Feynman-like approximation for the dynamic structure factor. This 
expression obviously belongs to the set of solutions of the (truncated) moment problem, i.e., 
it does satisfy all three frequency sum rules ([MS]). In a mathematical context, the function 
(fTTjl is usually referred to as a canonical solution of the moment problem. 

It must be pointed out that, unlike (fTTj) . the model expressions of Eq. (34) in [7] and 
Eq. (23) in [16] fail to satisfy the sum rules ([3]) and ([5|). In addition, for a strongly coupled 
system, F ^ 1, one might neglect the kinetic contribution, K (fc), in the moment C4 (fc). This 
leads to the inverse dielectric function obtained by Golden and Kalman by means of the QLC 
approximation [22]. 

As the Feynman-like approximation, expression (jlip reproduces a collective excitation at 
the frequency W2 (fc). The additional central peak here permits to account for hydrodynamic 
difussional effects [23] . Whereas in a multicomponent plasma this peak becomes dominant at 
low frequencies, in case of a one component plasma it vanishes as k decreases [24,25], as it can 
be observed from (jlip just by taking the long wave-length limit and recalling the well-known 
compressibility sum rule (see also expression P7)) below). In fact, this feature is related to 
the existence of a vanishing static (internal) conductivity in the long wave-length limit and 
might be understood as a macroscopic manifestation of the conservation of charged momentum 
during interparticle collisions at a microscopic level, which is obviously not applicable in a 
multicomponent system. 

Fast projectile polarizational stopping power. ~ Here we are interested in the fast pro- 
jectile stopping power of Coulomb systems with well-defined collective modes. Then it is 
sufficient to consider the canonical solution or the loss function in the form of (fTTjl . For a 
weakly coupled OCP, the coupling parameter, F or rg, is <C 1. Under such conditions one 
might make use of the RPA dispersion law which implies neglecting particle correlational 
contributions to the collective excitation frequency. Then, 



^2D 
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(12) 



Further, as a first approximation, one might consider only the asymptotic leading terms 
in the previous expression for the limits of large and small wave-vectors. Thus, 



■'2D 
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By substituting expressions (fTTjl and (fT^ into (fTj), one obtains 

^2 



dE 
dx 



VI - Sb^z- 



b^z^ 



dz. 



(13) 



(14) 



where we introduce the parameter 6 = {2kpaB) ^ and the new variable b = vp/v. The 
polynomial equation [16] z^ — b~^z + 6 = has one negative and two positive roots, Zmin, 
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fe« ,5-2/3. (15) 
The leading term of the asymptotic expansion of can then be calculated exactly: 



dx 



(16) 



This result was initially obtained by Bret and Deutsch [7]. 

On the other hand, if we want to consider all contributions in (|12p and the effects of an 
arbitrary coupling in the electron gas, i.e., to go beyond the RPA, we should consider tj| (fc) 
as given by pO|) . In our case, it suffices to consider the expansion of the expression (fTO)l into 
the power series in the wave-vector. For small k, we get [22] 



ul (fc i 0) ~ u. 



1 



Ef. being the correlation energy per particle. Similarly, for large k [22] 



w| {k t OO) ~ LU' 



2D 



^ I 3 , 3i?kin 
k + 



Thus, we might make use of the interpolation formula: 



k + 0{\) 



UJ2 (fc) ~ iiJ' 



1 



3£'kin , , 5Ec , 

Trk + — -k 



Again, by inserting expression (fTT|) . now with into formula ([1]) we get 

62 



dE _ 2V2Z^e^ 
dx rsttl 



where the term 



„ y/1 - - (562^-1 - 62^2 

, -Bkin , 5 Ec 



dz, 



(17) 



(18) 



(19) 



(20) 



(21) 



Ef 16 

accounts for the contribution of the average kinetic and correlation energies per electron. 
Clearly, this integral can be recast into 



dE 2V2Z2e2 



dx 



rsW 



=dz, 



(22) 



with 



&2 = 



1 - ^62 ■ 

Notice that in this case the condition (fTSl) transforms into: 



1 



(23) 



(24) 



The leading term of the asymptotic expansion of ((22|) for 6 — > obviously coincides with 
that of p6p . This remains true even if we abandon the singular form of the canonical solution 
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(jlip and substitute the (5-functions by adequate Gaussian densities, thus taking into account 
an exponentiaUy smaU damping of the plasmon mode. Indeed, the form (j22p is recovered then 
by a straightforward apphcation of the Laplace method of asymptotic analysis. 

Thus we prove that for the case of high- velocity projectiles, the correction to the stopping 
power due to correlational effects does not affect the main asymptotic term as v/v-p oo, 
obtained in Ref. [7]. Finite damping [22] can be taken into account using a non-canonical 
solution of the moment problem ([5]). 

Similar considerations are valid also for the 3D interacting electronic gas, where the Bethe- 
Bloch-Larkin [13-15] asymptote is not affected by correlations [26]. 

* * * 
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